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Abstract 

In energy harvesting (EH) relay networks, the coherent communication requires accurate estima¬ 
tion/tracking of the instantaneous channel state information (CSI) which consumes extra power. As a 
remedy, we propose two noncoherent EH relaying protocols based on the amplify-and-forward (AE) 
relaying, namely, power splitting noncoherent AE (PS-NcAE) and time switching noncoherent AE 
(TS-NcAE), which do not require any instantaneous CSI. We develop a noncoherent framework of 
simultaneous wireless information and power transfer (SWIPT), embracing PS-NcAE and TS-NcAE in 
a unified form. Eor arbitrary M-ary noncoherent frequency-shift keying (ESK) and differential phase- 
shift keying (DPSK), we derive maximum-likelihood detectors (MLDs) for PS-NcAE and TS-NcAE in 
a unified form, which involves integral evaluations yet serves as the optimum performance benchmark. 

To avoid expensive integral computations, we propose a closed-form detector using the Gauss-Legendre 
approximation, which achieves almost identical performance as the MED but at substantially lower 
complexity. These EH-based noncoherent detectors achieve full diversity in Rayleigh fading. Numerical 
results demonstrate that our proposed PS-NcAE and TS-NcAE may outperform the conventional grid- 
powered relay system under the same total power constraint. Various insights which are useful for the 
design of practical SWIPT relaying systems are obtained. Interestingly, PS-NcAE outperforms TS-NcAE 
in the single-relay case, whereas TS-NcAE outperforms PS-NcAE in the multi-relay case. 
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I. Introduction 

Conventional battery powered wireless eommunications systems require periodic recharging 
or replacement of the batteries, which incurs a high operation burden [HI and can be cumbersome 
or even impossible (e.g., for biomedical devices implanted in the human body |l2l). Recently, 
energy harvesting (EH) from the ambient radio-frequency (RF) signals has been developed as 
one of the attractive alternatives to prolong the lifetime of energy-constrained nodes in wireless 
networks [l3]|, flU. The dual use of RF signals for EH and information delivery (ID) has led to 
the novel architecture of simultaneous wireless information and power transfer (SWIFT) [|5l, ll^, 
which allows wireless nodes to scavenge energy as well as extract information simultaneously 
from the RF signals, thus constituting an appealing solution for energy-constrained applications 
such as wireless relay networks. 

Practical SWIPT receiver architectures make use of two different circuits performing EH and 
ID individually 0. The receiver may either switch between the EH and ID circuits in a time- 
division fashion, a scheme known as time switching |I8]|, or split the received RF signals into 
two streams fed to the EH and ID circuits at the same time, a scheme known as power splitting 
[|9l . Typically, EH and ID circuits operate with rather different receiver power sensitivities, 
i.e., —60 dBm for information receivers and —10 dBm for energy receivers ifTOl . [fTTI . Since 
the EH circuit is designed to maximize the EH efficiency while the ID circuit typically aims 
for maximum information rate, a fundamental rate-energy tradeoff exists for SWIPT systems 
0-1191, m- Various resource allocation and beamforming schemes were designed to achieve 
different tradeoffs between rate and energy in SWIPT systems. In particular, the optimum power 
allocation maximizing the information rate subject to EH constraint was studied for broadband 
systems with perfect instantaneous channel state information (CSI) llT3l . The optimum resource 
allocation with service differentiation between low- and high-priority data was studied for EH 
networks in [lT4l . In [lT5l . a joint optimization and game-theoretic framework was developed to 
optimize the packet delivery policy and cooperation strategy for delay tolerant networks with EH. 
Moreover, beamforming schemes which maximize the information rate subject to EH constraint 
or maximize the harvested energy subject to rate constraint were studied under imperfect CSI 
in [IT^ and lUTl . respectively. 

Recently, SWIPT has also been applied in wireless relay systems which allows energy- 
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constrained relay nodes to harvest energy from the souree RF signals and eonsequently brings in 
substantial benefits in wireless networks m. Various theoretieal analysis and praetieal design 
have been eondueted for EH relay systems in the literature. Speeifieally, for dual-hop amplify- 
and-forward (AF) networks, a greedy switehing protoeol between data relaying and EH was 
proposed in [fT9ll . and the joint relay seleetion and power alloeation seheme was developed in 
EOll . Moreover, the outage probability and ergodie eapaeity for EH AF relaying were derived for 
delay-limited and delay-tolerant applieations, respeetively, in [|2T]|. Besides AF relay systems, the 
EH has also been eonsidered for deeode-and-forward (DF) networks. In partieular, the outage 
probability of EH DE relaying was studied, taking into aeeount the spatial randomness of the 
souree-destination pairs, for various network topologies in Il22]l - [|25]| . 

The aforementioned studies on EH relay systems have either implieitly or explieitly assumed 
the instantaneous CSI availability to allow for coherent information deeoding. However, the 
eoherent SWIFT requires the souree to periodieally send training symbols, whieh ineurs an 
inereased signaling overhead and proeessing burden. Moreover, the relays may need to relay 
the training symbols and/or estimate the souree-relay ehannels ll26l . whieh results in additional 
energy eonsumptions and poses serious issues espeeially for energy-eonstrained relay nodes. As a 
remedy, the noncoherent SWIFT eliminating the need for the instantaneous CSI was first studied 
for the EH DE relay systems [|27]| . where two EH relaying protoeols, namely, the power splitting 
noneoherent DE (FS-NeDE) and time switehing noneoherent DE (TS-NeDE) were proposed, 
and the eorresponding maximum-likelihood deteetors (MEDs) faeilitating noneoherent SWIFT 
were obtained. However, these protoeols and the deteetors are applieable only for the DF relay 
systems. It is still unknown about how the noneoherent AF relaying ean benefit from EH, and 
how the noneoherent AF relaying performs as eompared to the noneoherent DF relaying in 
EH relay systems. In addition, for EH AF relay systems, the noneoherent MED minimizing 
the symbol-error rate (SER) in noneoherent SWIFT is still unknown. Realizing that none of 
the previous works have taekled the design ehallenges of noneoherent SWIFT in EH AF relay 
systems, we aim to fill the gap. 

In this paper, we study the noncoherent EH relay systems eonsisting of multiple AF relays, 
whieh ean harvest energy from the ambient souree signals and utilize the harvested energy to 
assist the eommunieation. The main eontributions of this paper are summarized as follows: 

• For arbitrary M-ary noneoherent frequeney-shift keying (FSK) or differential phase-shift 
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keying (DPSK) signalings, we propose a unified noncoherent SWIPT framework embraeing 
two EH relaying protoeols, namely, the power splitting noneoherent AF (PS-NeAF) and time 
switehing noneoherent AF (TS-NeAF). 

• Following the proposed SWIPT framework, we develop noneoherent MFDs for PS-NeAF 
and TS-NeAF in a unified form, whieh involves integral evaluations yet serves as the opti¬ 
mum performanee benehmark for noneoherent SWIPT. To avoid integral eomputations, we 
also develop a closed-form Gauss-Fegendre approximation based deteetor, whieh aehieves 
almost identieal SER to the MFD at substantially lower eomplexity. It is demonstrated 
that the proposed EH-based noneoherent deteetors aehieve full diversity in Rayleigh fading. 
In terms of the error performanee, the proposed PS-NeAF and TS-NeAF may outperform 
eonventional grid-powered noneoherent relay systems under the same total power eonstraint. 

• The ehoiee of the time switehing or power splitting parameters represents a tradeoff between 
EH and information transmission, and there exist unique optimum values of these parameters 
whieh minimize the SER. 

• The number of EH relay nodes is a key faetor on the performanee of the noneoherent EH 
relay systems: PS-NeAF outperforms TS-NeAF in the single-relay ease, whereas TS-NeAF 
outperforms PS-NeAF in the multi-relay ease. 

• The optimum relay position is eloser to the souree than to the destination, regardless of the 
EH relaying protoeols (This is in eontrast to the eonventional self-powered relay systems 
where the optimum relay position is eloser to the destination for AF and is eloser to the 
souree for DF). Furthermore, as the path-loss exponent inereases, the optimum relay position 
shifts slightly towards the destination, but is still eloser to the souree. 

• The M-FSK signaling with M > 8 is a more suitable solution for EH relay systems due to 
its higher energy effieieney eompared to M-DPSK, regardless of the EH relaying protoeols. 

• Comparison of the AF- and DF-based EH relaying protoeols demonstrates that the SER 
performanee of AF and DF relaying is almost the same in noneoherent EH relay systems, 
irrespeetive of the SWIPT arehiteeture (e.g., power splitting or time switehing). This is 
in eontrast to the eonventional self-powered relay systems where either AF or DF may 
outperform eaeh other depending on the relay loeations. 


The remainder of the paper is organized as follows. SeetionUUdeseribes the system model and 
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develops noncoherent EH relaying protocols. Section HII] obtains a unified noncoherent SWIFT 
framework for EH AE relay system. Section UV] derives the noncoherent detectors achieving 
SWIFT. Section |V] gains useful insights into the noncoherent SWIFT through simulations and 
Section |Vl] concludes the paper. 

Notation: We use (•)*, (•)^, (•)^, II ‘ II’ In(-) to denote the conjugate, transpose, 

conjugate transpose, expectation, real part, 2-norm, and natural logarithm, respectively. Also, 0, 
In, and in denote an all-zero column vector, an n x n identity matrix, and a column vector with 
1 at its n-th entry and 0 elsewhere, respectively. Moreover, A = B denotes that A is defined by 
B. Finally, x ~ CAf{fi, S) means that ai is a circularly symmetric complex Gaussian (CSCG) 
random vector with mean /x and covariance H. 

H. System Description and Noncoherent EH Relaying Frotocols 

Consider an EH relay network where the source terminal Tq communicates with the destination 
terminal T^ through the help of a set of potential relay candidates, where the selected relays can 
harvest energy from the source RE signals and forward the source information. Suppose that a 
total of K EH relays, T^, r = 1, 2, • • • ,K, are predetermined Q The source may communicate 
directly with the destination through the source-destination link hod or indirectly via the two- 
hop relay channels hor (first hop) and hrd (second hop). The relays, which do not have fixed 
power supplies, can only harvest energy from the RE signals radiated by the source and utilize 
that harvested energy to assist the source-destination communication through the noncoherent AF 
relaying, thus enabling SWIFT in the EH relay networks. We consider a composite fading model 
comprising the large-scale path loss Cij as well as the small-scale fading hij, ij G {Od, Or, 

The path loss component Cij is a distance-dependent constant and will be specified later. The 
small-scale fading coefficients are modeled as hij ~ CAA(0,1), which corresponds to the Rayleigh 
fading scenario. We assume that the instantaneous CSIs, hij, ij G {Od, Or, rd}^^, are unavailable 
to any terminals in the network, and thus, noncoherent signalings such as M-FSK and M-DFSK 
are considered. Note that the quadrature amplitude modulation (QAM) is inconsistent with the 


'Any potential relay selection protocols may be adopted. For example, a natural method is to choose the relays whose harvested 
energy is above certain threshold. Note that our proposed EH relaying protocols are valid for any specific relay selection schemes, 
and a detailed treatment of the relay selection scheme is beyond the scope of this paper. 
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Fig. 1. EH iC-relay systems with simultaneous wireless information and power transfer, where ID^ denotes the information 
delivery (or data relaying) from T^, fc = 0, • ■ • . jF. |(a)| PS-NcAF where each relay splits its received RF power into two portions: 
the p portion is for EH and the remaining 1 — p portion is for ID. |(b)| TS-NcAF where the total block time is divided into the 
EH phase of length aT and the ID phase of length (1 — a)T. 


noncoherent SWIPT framework eonsidered here as it typieally requires the instantaneous CSI at 
the receiver for eoherent information delivery. 


A. PS-NcAF Relaying Protocol 


In the PS-NcAF protoeol, eaeh EH relay node splits its received RF signal into two streams, 
whieh are fed to the EH and ID eircuits at the same time. Suppose that the total eommunication 
bloek time T (sec) is divided into K + 1 sub-blocks for PS-NcAF, each of length is alloeated 


to one of the transmitters (including the source and relays), as illustrated in Fig. |l(a)[ In the 
first sub-bloek, the souree broadcasts the RF signal with power Pq (Watts). At eaeh EH relay 
node, the p portion of the received RF signal is used for EH, and the harvested energy is used 
for relaying the souree signal in one of the following K sub-bloeks; thus, the average harvested 
power Pr available for data relaying at is 


Pr Tj pPqFqt ^ 


( 1 ) 


r = 1, • • • , Ji", where 0 < p < 1 is the power splitting faetor (PSF) and 0 < p < 1 is the EH 
effieieney ifT^ . The remaining 1 — p portion of the reeeived RF signal is fed to the ID circuit at 
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each relay node, where the signal is amplified using the harvested power in O and forwarded 
to the destination in one of the following K sub-blocks. 

1) PS-NcAF with M-DPSK: For M-DPSK transmission, each source message m e {0, • • • , M— 
1 } is differentially encoded into two consecutive information-bearing symbols s{l) and s(/ — 1) 
according to s{l) = s{l — where s(0) = 1 is the initial reference signal. Since 

the detection of each source message m is based on two consecutive received symbols at the 
destination, it is convenient to represent the signals as 2 x 1 vectors. Let s = [s(/ — l),s(/)]^ 
denote the source transmitted signal. Then, the received signal, y^r = [vori^ — 1), 2/or(0]^’ the 

relay node T^, r = 1, ■ ■ ■ , K, which is the intended for ID, can be expressed as 

yOr \/ (1 P'jPoTsP’Orhor^ “f \/l ~P'^0r A Vqj,, (2) 


where T* (sec) is the symbol time and the power-scaling factor 1 — p is due to the power splitting 
at Tr- The additive white Gaussian noise (AWGN) at (when Tq serves as the transmitter) is 
made up of two components: the AWGN due to the receive antenna (which is introduced before 
the power splitter), modeled as itor ~ CAA(0, ( 70 ^ 4 ^ 2 ), and the AWGN due to the ID circuit 
(which is introduced after the power splitter), modeled as ~ caa(o,< 2T2) a, m, El. 
At the destination, no power splitting is needed as the received signal is only used for ID. Thus, 
the received sinal at the destination can be written as 


Vod — 


\/ PoTgCodhodS + UQd + '^od, 


( 3 ) 


where uod ~ and v^d ~ CA/'(0, (yld, 2 l 2 ) are the AWGNs due to the receive 

antenna and ID circuit, respectively. 

Each relay amplifies its received signal yor with an amplifying gain Gr and forwards that 
signal to the destination. Thus, the received signal at the destination is given by 

yrd \/(l “ P)PoTsCQrCrdGr rd (a/1 - pUor + Vor) + Urd + Vrd, (4) 


where Urd ~ CA/’(0, al^ ^ 12 ) and Vrd ~ CAA(0, the AWGNs due to the receive antenna 

and ID circuit, respectively. The amplifying gain Gr ensures that the average transmission power 
for data relaying at T^ is fixed to Pr in ([T]) [|28]| . For PS-NcAF with M-DPSK, the amplifying 
gain Gr is given by 


Gr = 


P T 

1 rpj. g 


ypPoT^C 


Or 


E{bo.(/)P} 


(1 — p)PoTsCQr + (1 ~ P)^0r,l + G 


( 5 ) 


'Or ,2 
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2) PS-NcAF with M-FSK: For M-FSK transmission, the message m G {O,--- ,M — 1} 
is transmitted over one of the M orthogonal earners. The baseband equivalent reeeived signal 
intended for ID at the reeeiving terminal whieh is transmitted from Tj is denoted as an M x 1 
veetor = [yij(l), • • • ,yij(M)]^, ij G {Od, Or, Then, the signal model for PS-NeAF 

employing noneoherent M-FSK ean be represented as 


Yod — \/ PoTgCodhodSm+i + Uod + Vod, (6a) 

yOr \/(1 P^PoTgC-Q^hQ-r^rn+l T \/l "pUgr T Vgr, (6b) 

Yrd 's /(1 P)PoTsF'Or^rd(^rhorhrd^m+l T ^rd^rdi^y/^ PUQ,. T VQr) T T "Vrdi 

(6c) 


where Ujj ~ CA/'(0, and ~ CA/'(0, 2/m), ij ^ {Od, Or, rd})Lp are the AWGNs 

due to the reeeive antenna and the ID eireuit at Tj, respeetively. For PS-NeAF with M-FSK, 
the amplifying gain Gr, whieh ensures that the average transmission power for data relaying at 
Tr is fixed to Pr in ([U), is given by ll28ll 


Gr = 


P T 

± V <i 


r]pPoTsC 


■Or 


E{||yor|p} 


(1 — o)PnT.L, 








(7) 


B. TS-NcAF Relaying Protocol 


In the TS-NeAF protoeol, the total eommunieation bloek time T (see) is divided into two 
eonseeutive phases: the EH phase of length aT and the ID phase of length (1 —q;)T, as illustrated 
in Fig. |l(b)[ where 0 < a < 1 is the time switehing eoeffieient (TSC). In the EH phase, the 
souree sends RE energy flow to the relays with power Pq. Eaeh EH relay ean harvest energy 
from the reeeived RE signal and utilize that harvested energy for data relaying in the subsequent 


ID phase. The ID phase is further divided into K + 1 sub-bloeks of length 


K+l 


eaeh, where 


the first sub-bloek is assigned to the souree for new information transmission and eaeh of the 
remaining K sub-bloeks is assigned to one of the relays for data relaying. Then, the average 
harvested power Pr for data relaying at T^ is given by 

{K + l)r]PoCorO: 


Pr = 


a 


(8) 


r = 1, ■■■, K, where 0 < ?7 < 1 is the energy eonversion effieieney IfT^ . 
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1) TS-NcAF with M-DPSK: In the TS-NcAF protocol, the received signal during the ID phase 
at each node is solely used for information processing, i.e., no power splitting is involved. Thus, 
the baseband equivalent signal model for TS-NcAF with M-DPSK is represented in (2 x l)-vector 
form as follows: 


Vod — \/ PoTsCodhodS + UQd + VQdi 

Vor = y/ PoTgC-OrhorS + + Vq^, 

Urd 's/ PoTgC,Q^C,j.dGyhQj.hj.dS -\- -\- tlor) T '^rd T '^rd- 


(9a) 

(9b) 

(9c) 


r = 1 , • • • ,iT, where ~ CJ\f{0,afjil2) and Vij ~ CAA(0, 0 -^^^ 3 / 2 ), ij e {Od, Or, are 

the AWGNs due to the receive antenna and the ID circuit, respectively. For TS-NcAF with M- 
DPSK, the amplifying gain Gr, which ensures that the average transmission power at is fixed 
to Pr in dH]), is given by ll28]l 


Gr = 


{K + l)r]PoTsCora 


( 10 ) 


(1 — Q;)(PoT5£or + CTor.l + ^0r,2) 

2) TS-NcAF with M-FSK: The baseband equivalent signal model for TS-NcAF employing 
noncoherent M-FSK can be represented as 


Yod — y/ PoTsPodhodi-m+i + + vq^, (Ha) 

YOr 's/PoPsPor^oAm+l T ^Or T (Hb) 

Yrd PoTsPorPrdGrhQrhrd^m+1 T G^ \/Grdhrd (uor + Vor) + Urd + ^rd, (1 Ic) 


for r = 1, • • • ,K, where ~ CJ\f{0, and ~ CA/'(0, afj 2 lM), ij e {Od, Or, rdj^^, 

are the AWGNs due to the receive antenna and ID circuit at Tj , respectively. For TS-NcAF with 
M-FSK, the amplifying gain Gr, which ensures that the average transmission power at T^ is 
fixed to Pr in dH]), is given by 


Gr 


/ 


(1 — a 


{P + ^)vPoTsPor<^ 
) [PoTsPor + M{(^0r,l 


+ G 


Or, 2 


( 12 ) 


III. A Unified Noncoherent SWIPT Framework 


In the last section, we developed two EH relaying protocols, namely, PS-NcAF and TS-NcAF, 
which can be easily applied in conjunction with two widely adopted noncoherent signalings such 
as noncoherent M-FSK and M-DPSK, thus resulting in a total of four different noncoherent 
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SWIPT schemes for EH relay systems. From Section |nl we see that eaeh scheme has a different 
system model with a set of different parameters. In partieular, the amplifying gain expressions 
Gr eorresponding to the four possible sehemes in ([5]), ([7]), (fTOl) . and (fT2l) are all different. These 
distinet system models and inconsistent parameters make it rather eumbersome for further design 
and analysis. To resolve this problem, in this seetion, we will develop a unified SWIPT framework 
embracing both PS-NcAF and TS-NcAF, which enables unified further development. 

The key to the unifying proeess is to unify the definitions of all different system parameters. 
First of all, we introduee the effective noise variance for eaeh receiving terminal as follows: 


2 A 2 I 2 

^Od = ^Od,l + ^0d,2 

^2 A 2 I _2 

'^rd — ^rd,l “r ^rd,2 

2 A I ( 1 -P)c^ 0 r,l + 




(T, 


Or,l 


+ 


Or,2) 


PS-NcAF, 

TS-NeAF. 


(13a) 

(13b) 

(13c) 


Furthermore, let us introduce the information rate, R = ^ (bps), whieh is the total 

transmitted information bits normalized by the total eommunication bloek time T (sec), where 
Ng is the total number of transmitted information-bearing symbols. For PS-NeAF, the Ng new 
symbols are sent over the first sub-block of length -r^. For TS-NeAF, due to the time switehing 


effect, the Ng new symbols are transmitted over the first sub-block of length . By definition 
of R, we have 


log 2 ^ 


R = J {K+l)Ts ’ 

I (l-ajloggM 
I, iK+l)Ts ’ 

where Tg depends on T through Ng as follows: 


T - } iK+l)Ns-’ 
^ (1-«)T 


PS-NeAF, 

TS-NcAF. 


PS-NcAF, 

TS-NeAF, 


(14) 


{K+l)Ns 

For the purpose of performanee eomparison of PS-NeAF and TS-NcAF, one must ensure that 
the information rate R is the same for both sehemes, whieh is aeeomplished by ehoosing the 
eommon parameters T (total eommunieation bloek time), Ng (total number of symbols), and M 
(modulation alphabet size) for both sehemes. This implies that PS-NeAF and TS-NcAF must 
have different symbol duration Tg, in order to keep the same information rate R. 
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With the above definitions, the average signal-to-noise ratios (SNRs) of the direet souree- 


destination link, the souree-relay link, and relay-destination link ean be expressed, respeetively, 

ai] 


A 

lOd — 


A 

70r = 


A 


Pp^Od log2 ^ 
{K+l)((7‘^d,l+^Pd,2^^' 
{l-a)PoCod loft M 

(^+l)('^0d,l+'^0d,2)^’ 

{l-p)PoCor log2 M 
{K+l) [(l-p)o-g,. 
{l—a)PoCor log 2 AT 


PS-NcAF, 

TS-NeAF, 


J 


R 


PS-NcAF, 

TS-NeAF, 


ppPpCprCrd log 2 ^ 


\{l-p)PQCQr log2 M-|-(A' + 1) ((l-p)(Tg^^^-|-0-j 
{K + l)ar]PQCorCrd log2 ^ 


, 2 )«S] 7 


rd,l+<d,2) 


\{l-a)PQCor log2 J^+{K + l){cTl^^^+al^ .2)R^\ {^ld,l+^ld, 2 > ' 

for r = 1, • • • , iT, where .^ = 1 for M-DPSK and ^ = M for M-FSK. 


PS-NeAF, 

TS-NeAF, 


(15a) 


(15b) 


(15c) 


A. Unified PS/TS-NcAF Framework for M-DPSK 

Following the parameter definitions in (fT3]) - (fT5]) . the signal models for PS-NcAF and TS- 
NcAF employing M-DPSK can be expressed in a unified form as follows: 

VOr (^Or (\/ 'JOr^Or^ T "^Or); (16a) 

yod = o-Qd (^/^hodS + nod), (16b) 

yrd ^Qr^rdyJ^Or'^rdhQrhrd^r T ^Or^rd^J^rd^rdP^Or T ^rdP^rdi (16c) 

where riij ~ CN'{0, 12 ) for ij e {Od, Or, rd}f^^. Note that the unified model in (fT^ can represent 
either PS-NcAF or TS-NcAF by appropriately choosing the parameters according to (fT3])-(fT5]). 

B. Unified PS/TS-NcAF Framework for Noncoherent M-FSK 

Following the definitions in (fT3l) - (fT5l) . the signal models for PS-NcAF and TS-NcAF employ¬ 
ing noncoherent M-FSK can be unified as follows: 

yOr (POr (\/70r^0rlm-|-l T 11-Or) ; (I'^a) 


^The parameters 7 od and jor are defined as the actual link SNRs associated with the direct link and the first-hop link, 

respectively; but yrd is defined as the average SNR of the second-hop when the received signal from the first-hop is of unit 
A G^Crd 

energy, i.e., jrd = 7 ■ 
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(17b) 


yrd ^Or^rdy/'JOr'^rdhor “ 1 “ ^0r^rd'\/^rdhrdX^0r ^rdX^rdi 


(17c) 


where njj ~ CAA(0, Im) for ij G {Orf, Or, rd} and r = 1, • • • ,K. 

Remark 1: The main benefit of the unified noneoherent SWIFT frameworks in (fT^ and (fTTI) 
lies in the faet that it enables unified design and analysis for the PS-NeAF and TS-NeAF 
protoeols. For example, following these proposed unified noneoherent SWIPT frameworks, the 
noneoherent deteetors for both PS-NeAF and TS-NeAF ean be obtained in a unified form, whieh 
is treated in more detail in the next seetion. 


IV. Noncoherent Detectors for EH AF Networks 


In this seetion, the main objeetive is to develop the noneoherent deteetion sehemes for EH 
AE relay systems. To proeeed, we first taekle the mathematieal ehallenges involved in the 
development of the noneoherent deteetors. Then, the (exaet) MEDs for PS-NeAE and TS- 
NeAE are derived in a unified form, whieh eharaeterizes the optimum performanee benehmark 
for noneoherent SWIPT in EH AE networks. Einally, low-eomplexity noneoherent deteetors 
are derived, whieh aehieves almost identieal performanee to the MEDs at dramatieally lower 
eomplexity. 

A. Mathematical Preliminary 

The unified signal models in (fT^ and (fTV]) involve eomplieated transformations of the CSCG 
random variables/veetors. To faeilitate the development of the ME deteetion sehemes, it is useful 
to study the probability density funetions (PDEs) of those Gaussian transformations. 

Lemma 1: Consider random variables Xi ~ CJ\f{0, Hj), 2x1 random veetors Xi ~ CAA(0, crfl 2 ), 
and M X 1 random veetors yi ~ CAf{0, afiM), f = 1, 2, all of whieh are mutually independent. 
Eet 


ACq ~ X1X2C -f ^22:1 -f X2, 
To = XiX2ip -f X2yi + y2, 


(18) 


(19) 


where p is any integer number between 1 and M, c = [l,c]^, and c is any eomplex number. 



fxoix) 


( 20 ) 
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1 

[ ^ aj) 

' A ’ A 



1 


( 21 ) 


where x = [xi, X 2 Y, y = [ 2 / 1 , • • • , VmY, and /(ci, 62 , A, /^ 2 , A) is defined as 


p 

/(ci, 62 , A) A) A) = 


_ ( xA __I_ ^ 

+ 1 + 623:/ 


0 (1 + e2x)^(l + e2x ) 


dx, 


( 22 ) 


for any ei > 0, 62 > 0, A > 0) A > 0; A > 0. 

Proof: See Appendix ■ 

The PDF analysis in Lemma [T] is very useful for finding the MLDs for EH AF relay systems, 
whieh will be addressed in the next subsection. Furthermore, the generic analytical results 
obtained in Femma[I]may be useful for other different applications involving the same Gaussian 
transformations as considered here. 


B. Unified Noncoherent MLDs 


Following the unified noncoherent SWIPT frameworks in (fT^ and (fTTI) and the generic PDF 
analysis in FemmafU in this subsection, the MFDs for PS-NcAF and TS-NcAF are obtained in 
a unified form. 

Theorem 1: The MFDs for both PS-NcAF and TS-NcAF employing M-DPSK can be ex¬ 
pressed in a unified form as 


m = arg max 

1 1 -|- 2 yod 


27od ^{yoYi - 1)2/^,(/)+'2™/^} 


(Tn 


K / 

^ In J f al^'Jrd, (1 + 27or)(To^7rd, 

r> — 1 V 


Od J.—1 

IVrdil) - Vrdil - l)e^WM|2 _ l)ei2WM|2 






(23) 


Similarly, the MFDs for both PS-NcAF and TS-NcAF employing M-FSK are given in a unified 
form as 


j lOd |yod(’2r + l)P r/^ 2 , \ 2 

IlYrdlP - |yrd("i+ 1)P \yrd{m + l)\'^ 

Ad Ad 

Proof: See Appendix |Bl ■ 

The MFDs in (l23l) and (l24l) involve the computations of the integral /(ci, 62 , A; A, A) in 
(l22l) . which results in a high computational complexity. Nevertheless, the MFDs are very useful 


M- 1 


(24) 
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as they characterize the optimum performance benchmark for noncoherent SWIPT in EH AF 
networks. For example, one may use the MFDs as a theoretical performance upper bound to 
test the performance of any suboptimum detectors, if developed in the literature. 


C. Unified Noncoherent GLDs 


In this subsection, to avoid integral evaluations and make the MFDs suitable for practical 
implementation, we adopt the Gauss-Fegendre (GF) quadrature to approximate the integral, 
resulting in the so-called Gauss-Fegendre detectors (GFDs) in closed-form. 

Theorem 2: The GFDs for both PS-NcAF and TS-NcAF employing M-DPSK are given in the 
same form as (1^ . except that the multivariate function is replaced by 

which is given by 


J(ei, 62 , /3i, /32, A) — 


i=l 





(25) 


The parameters Wi and Zi are, respectively, the weights and nodes of the GF quadrature of order 
5, where wi = 128/225, W 2 = W 3 = m ;4 = ^5 = ^1 = 0, Z 2 = -Z 3 = 

2^5-2VW^, Z 4 ^ = -Z 5 = 2^5 -f 2^10/7, and ^{z) is defined as follows: 


f{z) = 


-( Pi : 132 ) 

0 \1 —ejlnz 1 — ) 

(1 — 61 In ^)^(1 — 62 In z) 


(26) 


Similarly, The GFDs for PS-NcAF and TS-NcAF employing M-FSK are given in (l2^ . where 
/(•, •, •, •, •) is replaced by /(•, •, •, •, •). 

Proof: See Appendix O ■ 

The MFDs in Theorem [2] require the computation of M integrals for detecting one information¬ 


bearing symbol, whereas the proposed GFDs avoid the integral computations and involve no 


special functions. Thus, the GFDs can be easily implemented in practice. Furthermore, we will 


demonstrate in the next section that the GFDs achieve almost the same error performance as the 
MFDs. Thus, the GFDs constitute practical solutions for EH noncoherent relaying systems. 


V. Numerical Results 

In this section, we evaluate the performance of the proposed noncoherent SWIPT schemes 
from various aspects through Monte Carlo simulations. 
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Fig. 2. SER versus PSF (p) or TSC (a) for PS-NcAF and TS-NcAF employing BDPSK and BFSK at information rate R = 1 
bps and SNR = 35 dB in a single-relay network (K = 1) with Dgr = 2. 

A A-relay line network model is eonsidered, where the relays are distributed over the straight 
line between the source and destination. Given the source-destination distance Dq^, the source- 
relay distance Dor and the relay-destination distance Drd may vary arbitrarily while guaranteeing 
Dor + Drd = Dod, fov r = 1, ■ ■ ■ , K. The EH efficiency is set to = 0.6 as in ifT^ . For ease 
of simulations, the antenna noise and ID circuit noise are assumed to have equal variances, 
i.e., afj i = 2 ij ^ {Od, Or, which ensures that each receiving node has a 

total noise variance of a^. Since the transmit power Pq applied to the source is the only energy 
supply to the whole network, the error performance of the whole system is parameterized by 
SNR = Pq/cq. Unless otherwise stated, the communication distance is set to Dod = 3 (m) and 
the bounded path-loss model ll2^ is adopted, i.e., Cij = ttW, where the path loss exponent is 
p = 4 as in [l25ll . We will evaluate the SER performance of the MEDs and GEDs for the proposed 
PS/TS-NcAF framework. In addition, the PS/TS-NcDF framework and the corresponding MEDs 
developed in IITtII will be considered as a benchmark for the purpose of comparison. 


A. Comparison of PS-NcAF and TS-NcAF 

In this subsection, we compare the error performance of PS-NcAF and TS-NcAF by investi¬ 
gating the impacts of the PSF p, the TSC a, and the number of relays K. The comparisons will 
be carried out for binary DPSK (BDPSK) and binary FSK (BFSK). 
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(a) 



(b) 

Fig. 3. SERs of PS/TS-NcAF and PS/TS-NcDF for multi-relay networks employing binary modulations. |(a)| 2-relay case (K — 2) 
with {Dor}r=i = {Ij 1-5}, R = 0.5 bps, and SNR = 30 dB. |(b)| 3-relay case (K = 3) with {Dor}r=i = {0.75,1.5, 2.25}, 
R = 1 hps, and SNR = 35 dB. 


1) Impacts of the PSF and TSC: The SER performance of the proposed EH relaying protocols 
is evaluated by varying the PSE p for PS-NcAE and the TSC a for TS-NcAE while fixing the 
SNR as SNR = 35 dB for binary (M = 2) noncoherent signalings at the information rate R = 1 
bps. Eig. [2] illustrates the SERs versus the PSE p or TSC a for the single-relay (K = 1) case 
where the relay is located at Dor = 1. We see that there exist (unique) optimum values of the 
PSE (i.e., p = 0.8) and the TSC (i.e., a = 0.4) which minimize the SERs of PS-NcAE and 
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Fig. 4. Grid-powered relay system with power Po/{K + 1) allocated to each terminal versus EFl relay system with Pq supplied 
to the source only for binary noncoherent signalings with R = 1 bps, where K = 3, {DoT-}r=i = {0.75, 1.5,2.25}, p = 0.8, 
and a — 0.55. 


TS-NcAF, respectively. This is because the choices of the PSF or TSC result in some tradeoff 
between EH and ID in EH relay systems. Specifically, for larger p (or a), more energy can be 
harvested and utilized for data relaying, whereas less signal energy (or time) can be used for 
ID. Eor smaller p (or a), more signal energy (or time) is available for ID; but less energy can 
be harvested and used for data relaying. Due to this tradeoff, the optimum values of p and a 
always lie between 0 and 1, which balances the operation of EH and ID such that the overall 
SER is minimized. 

2) Impact of the number of relays: Eor the single-relay network as illustrated in Eig. [2l the 
minimum SER of PS-NcAE (achieved around p = 0.8) is lower than the minimum SER of 
TS-NcAE (achieved around a = 0.4). However, this relationship is reversed in the multi-relay 
case. Specifically, for the 2-relay case (K = 2) employing binary noncoherent signalings where 
{ii)or}r=i = {1)1-5} and R = 0.5 bps, Eig. |3(a)| shows that the minimum SER of TS-NcAE 
achieved around a = 0.6 is slightly smaller than the minimum SER of PS-NcAE achieved 
around p = 0.85, for both BDPSK and BESK. In addition, for the 3-relay (K = 3) case where 
{ii)or}r=i = {0.75,1.5, 2.25} and R = 1 bps, as shown in Eig. |3(b)[ we can clearly see that the 


minimum SER of TS-NcAE achieved around a = 0.55 is much smaller than the minimum SER 
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of PS-NcAF achieved around p = 0.8. These comparisons indicate that the number of relays K 
is a key faetor whieh dietates to the superior EH protoeol. Speeifieally, PS-NeAF outperforms 
TS-NeAF in the single-relay ease {K = 1), while TS-NeAF outperforms PS-NeAF in the multi¬ 
relay case (K > 2). The reason for this is explained as follows. For PS-NcAF, since the EH 
time and the ID time are equal (both equal to the harvested power in ([T]) for data relaying 
at eaeh relay is independent of K, whieh makes the seeond-hop SNR 'jrd independent of K 
as well, as illustrated in (I15cl) . In eontrast, for TS-NeAF, due to the time switching operation, 
the EH time, aT, is generally unequal to the ID time, , at eaeh relay. Thus, for fixed 

harvested energy, the harvested power for data relaying at each relay becomes proportional to 
K + 1, and eonsequently, the seeond-hop SNR 'jrd is a monotonieally inereasing funetion of 
K, as illustrated in (I15el) . Henee, as K inereases, the seeond-hop link for TS-NeAF beeomes 
inereasingly more reliable as eompared to that for PS-NcAF, whieh results in extra performance 
improvement over PS-NeAF for iT > 2. 

For the multi-relay ease, we are also interesting in eomparing TS-NeAF and PS-NeAF at 
different SNRs, such that the performanee gain of TS-NeAF over PS-NcAF can be quantified at 
any target SERs. Fig. |4] shows the SERs versus the SNR for the 3-relay network employing the 
optimum EH parameters, i.e., p = 0.8 for PS-NeAF and a = 0.55 for TS-NeAF aeeording to 
Fig. |3(b)l where {Zi)or}r=i = {0.75,1.5,2.25} and R = 1 bps. For the target SER of 10“^, we 
see that TS-NeAF outperforms PS-NcAF by about 1.5 dB, for both BDPSK and BFSK. This 
signifieant gain is very attraetive as the energy effieieney is of major eoneern in EH systems. 
Furthermore, the eonventional grid-powered relay system is evaluated as a benehmark. For a fair 
comparison in terms of the total power consumption, the powers allocated to the souree and all 
relays are equal to ensuring the same total power of Pq as for the EH relay system. We 
observe that both EH protoeols, namely PS-NcAF and TS-NeAF, outperform the grid-powered 
AF relaying protocol, for the same noncoherent modulation. This is because given the same total 
power, the EH ability allows the relays to obtain extra powers from the ambient RF signals, thus 
boosting the overall performanee over the grid-powered relay system with no EH. 


B. Impacts of Relay Positions and the Path Loss Exponent 

1) Impact of the relay position For eonventional self-powered relay systems where the 
EH eapability is disabled l[29]l - lf34]| . the impaet of the relay position on the error performanee 
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(b) 

Fig. 5. SER versus the relative relay position {Dor/Dod) for PS/TS-NcAF and PS/TS-NcDF in the single-relay case with PSF 
p = 0.8, TSC a = 0.4, and information rate R = 1 bps. |(a)| g = 2.7 and SNR = 26 dB. |(b)| p = 4 and SNR = 35 dB. 


is well understood. Specifieally, the optimum relay position minimizing the SER in coherent 
relay systems is at the midpoint between the souree and destination, for both AF and DF [|29ll . 
If^ . For noncoherent relay systems, however, the optimum relay positions differ for AF and 
DF. For DF relaying, the optimum relay position is eloser to the souree than to the destination 
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on . on . whereas the optimum relay position for AF relaying is eloser to the destination than 
to the souree On . 04ll . 


Fig. [5] illustrates the SER versus the relative relay position Dq ^/for noneoherent EH relay 
systems. For both AF-based PS/TS-NeAF proposed in this paper and the DF-based PS/TS-NeDF 
in IIITII, the optimum relay positions are Dor = 0.2Dod when the path loss exponent is g = 2.7, 


as shown in Fig. |5(a)[ and D^r = O.SDod when the path loss exponent is ^ = 4, as shown 


in Fig. |5(b)[ That is, unlike the eonventional noneoherent relay networks where the optimum 
relay positions depend on the relaying protoeols (AF or DF), the optimum relay positions for 
noneoherent EH relay systems are invariant with respeet to AE and DE, and are always eloser 
to the souree than to the destination. This major differenee is due to the faet that EH relays are 
solely powered by the source, and thus, the dominant performance limiting factor in EH relay 
networks is the harvested energy at the relays, which is invariant with respect to AE and DE. 
Moreover, the overall performance is affected by the path loss, which is a distance-dependent 
constant. Since both the EH operation and the path loss effect are irrelevant to the relaying 
protocols, the optimum relay positions are the same with respect to AE and DE in EH relay 
systems. We note that the EH operation and the path loss have different impacts on the relay 
position. On one hand, the relays must be sufficiently close to the source in order to harvest 
enough energy. On the other hand, the relays must not be infinitely close to the source (e.g., 
located at the source) because otherwise the relay-to-destination path loss would become very 
significant, which degrades the overall performance. Considering both the EH and path-loss 
effects, the optimum relay position for EH relay systems is relatively closer to the source than 
to the destination. 


2) Impact of the path loss exponent g: As illustrated in Eig. [51 the optimum relay position, 
which is always closer to the source than to the destination, shifts slightly towards the destination 
as the path loss exponent g increases from 2.7 to 4. This is because a larger path loss exponent 
corresponds to a larger relay-to-destination path loss, which can degrade the overall error perfor¬ 
mance. As such, the relay should step towards the destination a little bit to diminish the effect 
of the increased path loss. Nevertheless, the prerequisite of EH operation cannot be obviated, 
and thus, the relays must always be deployed closer to the source than to the destination, with 
some slight movements towards the destination as the path loss exponent increases. 
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(a) 



(b) 

Fig. 6. SERs of PS/TS-NcAF and PS/TS-NcDF for higher-order constellations in the single-relay case {K = 1) with Dor = 1. 
[(a)|M = 4, SNR = 38 dB, and i? = 2 bps.j^M = 8, SNR = 40 dB, and R = 3 bps. 


C. Impact of Higher Constellations 


So far, we have only eonsidered noneoherent binary signalings with M = 2, and a 3 dB 
performanee gain of BDPSK over BFSK ean be observed in Fig. |4l In this subseetion, we 
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Fig. 7. SER versus SNR for PS-NcAF and PS-NcDF employing BDPSK and BFSK: 1) single-relay (K — 1) case with 
Dor = 2, p = 0.8, and R — 1 bps and 2) two-relay {K = 2) case with {Dor}r=i = {li 1-5}, p = 0.85, and R = 0.5 bps. 


will evaluate the performanee of higher order noneoherent signalings for M > 2. Speeifieally, 
for a single-relay network with Dq^. = 1, the SER of M = 4 at SNR = 38 dB and R = 2 


bps is illustrated in Fig. 16(a)i where 4-DPSK slightly outperforms 4-FSK, irrespeetive of the 
EH relaying protoeols (e.g., PS/TS-NeAF or PS/TS-NeDF). Moreover, the SER for M = 8 at 
SNR = 40 dB and R = 3 bps is illustrated in Fig. 16(b)l which shows that 8-FSK substantially 
outperforms 8-DPSK, irrespective of the EH relaying protocols. Observing the trend from M = 2 
to 4 and 8, we conclude that M-FSK is more efficient in terms of energy consumption than M- 
DPSK, which results in significant performance gain for M > 8. In particular, the performance 
gain of 8-FSK over 8-DPSK is about 4 dB in the target SER of 10“^, as demonstrated in 


D. Comparison of PS/TS-NcAF with PS/TS-NcDF / l27l/ 

Since both AF and DF relaying protocols are useful candidates for achieving SWIPT in 
EH relay systems, it is meaningful to compare the performance of the AF-based PS/TS-NcAF 
proposed in this paper with that of the DF-based PS/TS-NcDF developed in ll27]l . The SER 
performance of the EH relaying protocols generally depend on the relay positions, the EH 
parameters (such as PSF and TSC), and SNRs. Thus, performance comparisons can be carried 
out from various angles. In particular. Fig. [5] illustrates the comparison of PS/TS-NcAF and 
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Fig. 8. SER versus SNR for TS-NcAF and TS-NcDF employing BDPSK and BFSK: 1) single-relay (K = 1) case with 
Ddr = 2, a — 0.4, and R = 1 bps and 2) two-relay {K = 2) case with {Dor}r=i = {L 1-5}, Q = 0.6, and R — 0.5 bps. 


PS/TS-NcDF by varying the relay positions for fixed EH parameters and fixed SNR. In addition, 
Fig. [^illustrates the SERs for a variety of EH parameter settings (PSF or TSC), where the relay 
position and SNR are fixed. The remaining possible comparison is therefore to fix the relay 
position and EH parameters, while varying the SNRs. To this end, we compare PS-NcAF with 
PS-NcDF, and TS-NcAF with TS-NcDF at different SNRs in Figs. |7] and [8l respectively, where 
two typical network geometries are considered: 1) single-relay (K = 1) case with Dqj. = 2, 
a = 0.4, p = 0.8, and R = I bps and 2) two-relay (K = 2) case with {Dor}r=i = {1,1-5}, 
a = 0.6, p = 0.85, and R = 0.5 bps. Note that the PSF p and TSC a employed in the comparisons 
of Figs. |7] and [8] are optimum values determined in Figs. [2] and |3(a)[ For all the comparisons 
carried out from various angles in Figs. [5]-[8l we observe that the SER curves associated with 
PS-NcAF and PS-NcDF almost overlap, so do the SER curves associated with TS-NcAF and 
TS-NcDF. These observations reveal that the noncoherent AF and DF relaying protocols yield 
almost the same error performance in EH relay networks, regardless of the SWIPT architectures 
(e.g., power splitting and time switching), the modulation types (e.g., DPSK and FSK), and 
system parameters (e.g., information rate R, number of relays K, relay positions, and SNRs). 
This is in contrast to the conventional self-powered relay networks, where either AF or DF may 
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outperform each other depending on the relay positions. In particular, for self-powered relay 
networks, DF might outperform AF when the relay is much closer to the source than to the 
destination, because the noise amplification due to AF substantially degrades the performance. 
On the other hand, AF might outperform DF when the relay is much closer to the destination, due 
to the high relay decoding error probabilities in DF. However, both effects of noise amplification 
(associated with AF) and error propagation (associated with DF) are significantly weakened and 
become less important in EH relay systems, because the dominant performance limiting factor 
is the harvested energy at the relay, which is invariant with respect to the AF and DF relaying 
protocols. Hence, for EH relay systems, there is almost no difference between the noncoherent 
AE and DE relaying in terms of the error performance. 

In Eigs.|7]and[8l we also observe that the proposed noncoherent detectors achieve full diversity 
orders of 2 and 3 for the single- and two-relay networks, respectively, with the direct link. 
Eurthermore, we would like to note that the SER performance of the GLD is in excellent 
agreement with the SER performance of the MED, as demonstrated in all Eigs. [2]-l8l Moreover, 
the GLD is given in closed-form, involving no integrals or special functions at all. This suggests 
that the closed-form GLD is an appealing practical solution for noncoherent SWIFT in EH 
AE relay systems. Einally, this also indicates that J(ei, € 2 ,/9i,/ 92 , A) of (1251) is a very accurate 
approximation for the integral J(ei, 62 , /Si, /S 2 , A) in (l22l) . 

E. Performance Evaluation for Practical Indoor Communications Scenario 

We now evaluate the performance of the proposed detectors under a practical indoor communi¬ 
cation scenario operating at 900 MHz, where the typical communication distance of Dod = 10 (m) 
is considered. Note that according to the current state-of-the-art of RE energy harvesting/transfer, 
the maximum line-of-sight (LOS) operating distance for some RE energy harvesting chips is 12- 
14 meters [l35ll . If3^ . In this paper, we consider Rayleigh fading with blocked LOS due to a 
double plasterboard wall; thus, the communication distance of 10 meters is a very reasonable 
and practical choice. According to empirical measurements carried out at 900 MHz,, we set 
the path loss exponent as ^ = 1 . 6 , which corresponds to the communication inside the typical 
office building ll37l Table 2.2]. The partition loss is set to 3.4 dB assuming that the transmit and 
receive terminals are separated by the double plasterboard wall ll37l Table 2.1]. Also, we assume 
that the communication occurs on the same floor; thus, the floor loss can be ignored. According 
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Fig. 9. SER performance for typical indoor communications with path loss Lij (dB) = lOlogj^g L ^ —3.4, where D^d ~ 10 
(m), Dor = O.lDor for r = 1, 2, 3, a = 0.55, p = 0.8, and R — 1 bps. 


to the indoor path loss model featuring at 900 MHz 11371 eq. (2.38)], we set the path loss as 


Cij (dB) = 10 logi 


-ij ^^610 i+Di fc> 

' ^3 


3.4 in this subsection. Fig. [9] shows the SER performance for an EH 


relay system with K = 3 relays, all located at Dq^. = OADq^, r = 1, 2, 3. Note that we choose 
this relay position because the optimum relay position is much closer to the source than to the 
destination, as indicated in Eig. [31 Also, the EH parameters are chosen as a = 0.55 and p = 0.8, 
which are optimum values for K = 3 relays according to Eig. |3(b)[ Our simulation in Eig. 
demonstrates that for binary modulations with R = 1 bps, the diversity order of 4 is achieved 
for iF = 3 EH relays (with the direct link), which confirms the full diversity performance. 


VI. Conclusions 

In this paper, we have developed a noncoherent SWIPT framework for EH AE relay systems, 
which embraces both PS-NcAE and TS-NcAE in a unified form and supports arbitrary M-ary 
noncoherent signalings including M-ESK and M-DPSK. The main advantage of the proposed 
SWIPT scheme is that it eliminates the need for the instantaneous CSI, which alleviates the 
system overhead and reduces the energy consumption. Eollowing this framework, we developed 
noncoherent MEDs for PS-NcAE and TS-NcAE in a unified form, which involved integral 
evaluations yet served as the optimum performance benchmark for noncoherent SWIPT. To avoid 
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integral computations, we also proposed closed-form GLDs, which achieved almost identical SER 
performance to the MLDs at substantially lower computational complexity. It was demonstrated 
that these noncoherent detectors achieve full spatial diversity in Rayleigh fading. Moreover, in 
terms of the error performance, the proposed SWIFT relaying schemes may outperform con¬ 
ventional grid-powered relay systems under the same total power constraint. Finally, numerical 
results led to useful design insights into the noncoherent SWIFT in various aspects, including 
the effects of the time switching or power splitting parameters, relaying protocols, the number 
of relays, relay positions, and the modulation alphabet size, which are summarized as follows: 

• The choice of the power slitting factor or time switching coefficient represents some tradeoff 
between energy harvesting and information delivery. Unique optimal values of the FSF or 
TSC minimizing the SER exist between 0 and 1. 

• The number of EH relays is a key factor on the performance of noncoherent SWIFT: 
FS-NcAF outperforms TS-NcAF in the single-relay case, whereas TS-NcAF outperforms 
FS-NcAF in the multi-relay case. 

• The optimum relay position is closer to the source than to the destination, with some slight 
movements towards the destination as the path loss exponent increases. 

• The noncoherent M-FSK with M > 8 is more energy-efficient and thus might be more 
suitable than M-DFSK for EH relay systems, regardless of the EH relaying protocols. 

• The noncoherent AF and DF relaying protocols yield almost the same error performance 
in EH relay networks, regardless of the SWIFT architectures, the modulation types, and 
system parameters. 

An interesting extension of this work is to analyze the SER of the MFD or GFD, and using 
the obtained SER expression to analytically determine the optimum time switching coefficient or 
power splitting factor. In addition, more sophisticated and intelligent coordination of the source 
and relays’ transmissions (such as relay selection) may also be considered in order to further 
boost the performance. 


Appendix A 
Froof of Femma[I] 

In this proof, we will first study the distribution of Xq, followed by a brief derivation of the 


FDF of Yo. 
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It is not hard to show that the conditional distribution of Xq given IX 2 P = 2 : is Xo||X 2 | 2 = 
CA/'(0, So), where 

(f 2 i + (jf)z + cr|, Qic*z 


Sn = 


Qicz, 


(fli|cp + af)z + (T2 


(A.l) 


We denote by det(So) the determinant of Sq, and it can be shown that 

det(So) = (cr^z + (T 2 )([f^i(l + |cp) + a^]z + 

1 


So-' = 


det(So) 


{D.i\c\'^ + af)z + a2, —ilic*z 

—VLicz^ (f2i + < 7 ^) 2 ; + cr| 


(A.2) 

(A.3) 


It follows that the conditional PDF of Xq given IX 2 P = 2 : is / 


Xo\\X 2 \^=z^^^ 7r2det(So) 


Substituting in the expressions of (IA.2I) and (IA.3I) . we can show that 




tC q tC — , ^ ^ ^ *4- 2 


* * 1 
X^,Xn 


(f2i|cp + ( 7 ^) 2 : + (72, —f2iC*2: 

—D1C2;, (Di + a‘l)z + 


det(So) 

[(Di|cp + 0 - 1 ) 2 ; + o-|]|a;ip + [(Di + 0 - 1 ) 2 ; + o-2]|a;2p - 2Di2;5R(ca;ia;2) 


Xi 

X2 


Ai 


+ 


{alz + o-|)([Di(l + |c| 2 ) + afjz + al) 
A 2 


0 -^ 2 ; + o-| [Di(l + |cP) + o-f]2; + o-|’ 

where (IA.6I) follows by the partial fraction expansion. It can be shown that 

\X2 - CXip 


Ai = ( 0 - 12 ; + o-2)a;^So 


-2 “ 


(A.4) 

(A.5) 
(A.6) 

(A.7) 
(A.8) 

(A.9) 

Jo (0-^2; + 0 -|)([fli(l + |cP) + 0-1)2; + 0-2) ^ ^ 

Applying the change of variable ^ t in (lA.lOl) and taking some algebraic manipulations, we 
can show that 


1 + |cP 

“ 1 

A 2 = ([^ 1(1 + IcH + 0 - 1 ) 2 ; + o- 2 )ic'^So 


\xi + C*X2p 
1 + |c|' 


ni(i+|c|2)+<T^ 

The unconditional PDF of Xq can be computed as follows: 

/•oo 


fxoix) = / / 


,ix)f\x^\2{z)dz 


'0 


Xo||X2|2 = 

-(^+ -^1 
2 ^00 ^ \Q2 


2 + (T^ ' [n]^(l + |c|2) + (T^]2 + (T^ 


) 


-dz. 


exp 


fxo{x) = 


— f + 




B2 


n,aT“ + . n2[»i(i+kP)+'^il 


1+^^i 1+ 


(7rcT|)2 Jq ^ ^ 


-d2, 


(A.ll) 
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where 


A ^1 _ \X2- CXip 

' ~ ~ (1 + IcHcri 


^2 


2 ’ 


(A.12) 


\Xi + C*X2\ 


al (l + |c|>r 

Finally, following the definition of /(ci, 62 ,/9i,/ 32 , A) in (l22l) . we ean simplify (lA.llI) into (l20l) . 

The derivation of the PDF of Yq follows similar lines as the derivation of /xo(*)> whieh is 
summarized in two main steps. In the first step, the eonditional PDF of Yq given IX 2 P = 2 : is 
obtained as 


\yp\' 


M 


Ij/iC 


g (fll+a'flz + a':^ g a'fz + a'^ 

•^^011X212=2^^^ 7 r[(f 2 i + a^z + cr^] vr(cr^ 2 : + cr|) 

i¥=P 


\ypr 


(Qj^+o-^)2+cr^ 




7r^((T^2; + cr|)^ ^[{0.1 +al)z + al] 

In the second step, the unconditional PDF of Yq is obtained as 


poo 


exp 


- u + 


\?-\yv\^ I 
ZT + 


Q.2crfz+a2 


\yp? 

fl2{fZli+cr'f}z+cr2 


(TTcxi) 


2 \M 


l + ^z 


M -1 r 


1 I 02 (Oi+o-?) „ 
J- -I- „2 2 


dz, 


which simplifies to (1211) following the definition of /(ei, 62 ,/Si,/ 92 , A) in (l22l) . 


(A.14) 


(A.15) 


Appendix B 
Proof of Theorem [I] 

In this proof, we will start with the derivation of the MLD for M-DPSK. Then the MLD for 
M-FSK will be derived in a similar procedure. 

Recall the unified signal model in (fT^ for both PS-NcAF and TS-NcAF employing M-DPSK. 
It is not hard to see that yod\m ~ CJ\f{0, Sod), where 


Sod — <Xi 


Od 


1 + 7od 
lod^ 

Thus, the conditional PDF of yod given m can be expressed as 


—2‘1'KvnlM 


7ode 

J2™/M ^ ^ 


(B.l) 


M exp I , I, (3,,, 


+ 27od) 


1 + 27o(i 


cr, 


Od 


/■ 
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Using some algebraic manipulations, the signal model for Urd in (libel) can be rewritten in the 
desired form as Xq in (fTSl) 


yrd = X 1 X 2 C + X 2 X 1 + X 2 , 


(B.3) 


where 


c=[i,c{i)r 

(B.4a) 

Xl = aQr^/^hors{l - 1) ~ CA/'(0, cro^7or) 

(B.4b) 

X 2 ^rd\/^rdhrd ~ CA/"(0, ^rdH^d^ 

(B.4c) 

Xl = aor-nur ~ CAA(0, al^h) 

(B.4d) 

X 2 = (Ird-rird ~ CJ\f{0, (Jrdh)- 

(B.4e) 


Applying Lemma [H the conditional PDF of yrd given m can be written as 




fiyrd\m) = \ I ( a^r'yrd, (1 + 27or)cTo^7rd, 


\yrd{l)-yrd{l-l)ei^^^^^\‘' 




\yrd{l) + yrd{l-l)ei^^^^^? 


2^rd 


,1 • 


Following the definition of the MLD, we have 


m = arg max f{yod, {yrd}r=i\m) 

m= 0 ,-" ,M —1 


K 


= arg max In fiyodlm) + \nf{y 

t r=l 


’rd\m) 


(B.5) 


(B.6) 

(B.7) 


where f (yod, {yrd}^=i\'i^) is the likelihood function when m is transmitted, and (IB.7I) follows 
by the independence among all different signal branches. Substituting (IB.2I) and (IB.5I) into (IB.7I) 
yields the MLD in (1231) for M-DPSK. 

The MLD for M-FSK can be derived with similar steps as in the derivation of the MLD for 
M-DPSK. Specifically, we first obtain the PDF of yod in (I17bl) . conditioned on the transmitted 
message m, as 


f(yodjm) = 


llyodir 

e ""od 




T'od |yod("i+i)l 

, l+'TOd ^ 


(B.8) 


Then, yrd in (I17cl) can be rewritten in the desired form as l^o ns yrd = XiX 2 im+i + -^ 2^1 + y 2 , 
where Xi = aory^hor ~ CJV(0, a^^Jor), X 2 = (^rdy/^Kd ~ CN'{0,af^yrd), Vl = (^Or^Or ~ 
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CAf{0,a^jM), and ^2 = (^rd^rd ~ CAf{0,a^jM)- Then, applying Lemma [H the conditional 
PDF of Yrd given m can be expressed as 



Finally, following the definition of MLD, we can show that the MLD for M-FSK is given by 

dal. 


Appendix C 
Proof of Theorem [2] 


The main difficulty involved in the MLD is the integral computation of J(ei, 62 ,/^i,/ 92 , A) in 
(l22l) . Thus, we will develop a tight closed-form approximation for this integral. To that end, we 
first transform the original integral into a desired form which is easy to approximate. Specifically, 
we carry out the following algebraic operations: 



/(ci, 62 , /3i, {32, A) 


(C.l) 


(C.2) 


(C.3) 


where (1C.21) follows by the change of variable z and ^jJ{z) is defined in (|2^ . The 

integrand ^jJ{z) is a well-behaved function: first, it is a bounded positive function with 0 < 
^jJ{z) < 1 ; second, it is a smooth function with continuous derivatives with respect to 2 ; of any 
order. Furthermore, the integral limits are finite. These elegant properties enable very accurate 
approximation of the integral using the GL quadrature with a finite order. Indeed, the choice of 
the approximation order must always compromise between the accuracy and computational cost 
lf38]l . The higher order, the more accurate the GL quadrature is. On the other hand, however, 
higher order requires the computation of more terms which involves higher computational cost. 
Considering both the accuracy and computational cost, we choose the fifth order GL quadrature 
for our purpose (Ml, which results in J(ei, 62 , 13i, (32, A) /(ei, 62 , /Si, /S 2 , A), where /(•, •, •, •, •) 
is given in closed-form in (l25l) . Finally, using in place of in the MLDs 
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yields the corresponding GLDs 
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